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NOVEL DIFFUSIVITY MEASUREMENT TECHNIQUE 


Nasser Rashidnia 

National Center for Microgravity Research 
Cleveland, Ohio 44135 


SUMMARY 

A common-path interferometer (CPI) system was developed to measure the diffusivity of liquid pairs. The CPI 
is an optical technique that can be used to measure changes in the gradient of the refraction index of transparent 
materials. This system uses a shearing interferometer that shares the same optical path from a laser light source to 
the final imaging plane. Hence, the molecular diffusion coefficient of liquids can be determined using the physical 
relations between changes in the optical path length and the liquid phase properties. The data obtained with this 
interferometer were compared with similar results from other techniques and demonstrated that the instrument is 
superior in measuring the diffusivity of miscible liquids while keeping the system very compact and robust. CPI can 
also be used for studies in interface dynamics and other diffusion-dominated-process applications. 


INTRODUCTION 

Optical diagnostic techniques have become an integral part of many measurement applications in industrial and 
research laboratories. In recent years, the point diffraction interferometer (PDI) and its phase-shifted version have 
been used to measure optical wave fronts for lens testing and combustion and fluid flow diagnostics (refs. 1 to 3). 
The PDI is considered robust (ref. 4) because it has a common-path design. Another technique, the schlieren, has 
been used for many years for flow visualization and gives an overall qualitative picture of the variations in the index 
of refraction within a transparent test section. However, there are advantages in using the shearing interferometer, 
which has the same capability as the schlieren head. The CPI allows one to easily align the focused light on the 
interferometer head, and it can be used instead of a schlieren head and a PDI to make quantitative measurements in 
liquids that have large variations in the index of refraction, as is often the case in interface dynamics studies. This 
CPI capability is due to its basically being a Wollaston prism in combination with a polarizer and an analyzer; thus, 
to measure large variations in the index of refraction, one need only to use different Wollaston prisms for each par- 
ticular experimental condition. 

The dynamics of miscible interfaces is an active area of research that has been identified to benefit from experi- 
mentation in reduced gravity. Recognized as an important transport phenomenon is the study of the flow patterns by 
the interface when one liquid is slowly displaced by another such that diffusion plays an important role in the 
dynamics. It has been suggested that nontraditional stresses in the fluids, resulting from the steep variation in the 
concentration of the miscible fluids in the mixing zone, might be important in the dynamics. These effects are over- 
whelmed by the buoyancy-caused flow under terrestrial conditions. The diffusion coefficient (diffusivity) of the 
miscible fluids is a property that is important to the proper design of experiments. In particular, it determines the 
desirable range for the speeds of the displacing fluids so that the effect of diffusion is not overshadowed by convec- 
tive transport of the mass. Petitjeans and Maxworthy (ref. 5) used a variation of Wiener’s technique, described in 
Sommerfeld’s textbook on optics (ref. 6), to measure the diffusion coefficient. A narrow laser beam was passed 
through a quartz cell in which two liquids were diffusing into each other across a horizontal miscible interface. As 
the beam traversed the test cell, it was bent vertically when it traversed a vertical gradient of concentration. A light 
detector was then used to measure the deflection of the exiting beam as a function of vertical position. The deflec- 
tion angle is related to the vertical gradient of the index of refraction and hence to the concentration gradient within 
the test cell. The average diffusion coefficient was obtained from an analysis of the measured concentration gradient 
profile. The analysis assumed that the diffusion process is one dimensional and is characterized by a constant value 
of the diffusion coefficient. The technique used by Petitjeans and Maxworthy (ref. 5) to measure the concentration 
gradient is a point measurement that requires the light beam and the detector to be traversed through the region of 
interest. Therefore, instantaneous measurement of the concentration gradient cannot be obtained everywhere near 
the interface. This limitation can be overcome by using a CPI technique that employs Wollaston prisms to obtain the 
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distribution of the index of refraction in the mixing region. The specific technique incorporates a Wollaston prism 
operating in finite fringe mode (ref. 7) to visualize the concentration gradient in the mixing zone in real time. This 
measurement requires a novel approach in setting the prism at a particular angle with respect to the mixing interface 
orientation. The diffusivity is then obtained from the measured concentration profiles. The author used the finite 
mode of fringe arrangement of the CPI (Petitjeans and Maxworthy’s technique, ref. 5) and reports herein the meas- 
ured diffusivity of three combinations of pure silicone oils. 


DESCRIPTION OF CONCEPT AND MEASUREMENT APPARATUS 


Interferometric Model Concept 

The optical setup of many shearing interferometers, including the subject CPI, is similar to that of a schlieren 
system; therefore, the term “schlieren interferometer” has become popular. Among a great number of interferometer 
systems, an arrangement with a Wollaston prism as the shearing element is most frequently used (ref. 7) in experi- 
mental fluid mechanics, gas dynamics, and aerodynamics. This discussion will concern the customized CPI system 
used for the diffusion coefficient measurements reported herein. In contrast to the schlieren technique, the signal 
obtained with a CPI results from the change in the optical phase of light waves passing through the transparent 
object, not by light deflection. Two parallel, laterally sheared light rays passing through the test section at coordi- 
nates (y + d/2) and (y - d/2) are considered (ref. 7), where d is a small distance between the two rays. The shearing 
element positioned in or near the focal point of the schlieren head ensures that the two rays coincide after they have 
passed through this element. Under the assumption that the conditions of optical coherence are fulfilled, the two 
coinciding light rays can interfere with each other. If the refractive indices n (or liquid concentration c ) at coordi- 
nates (y + d/2) and (y - d/2) in the test section are not equal, two light rays (or waves) exhibit a difference in their 
optical paths Al(xy) as given by 


Z2 

A l(x,y)= | n(x,y + — 9 z)dz — 

Zl 


J r^x,y-^,z^dz 


( 1 ) 


Two assumptions have been made for deriving this equation: first, the coordinates where the rays enter and leave the 
test section z j and z 2 are equal for the two rays; that is, the test section is bounded by plane walls. Second and more 
important, the propagation of light in the test section is straight and not deflected. This assumption is reasonable for 
most gas flows and for very low to moderately sheared liquids. One should be aware that this assumption might be 
violated in some liquid flows and gases with very strong density gradients (e.g., in plasma flows). Exceptions from 
the validity of this assumption, known as “strong refraction,” will not be of concern here. Strong refraction is much 
more complicated because neither the refractive index (the concentration or density) nor the path of the light in the 
test section is known. 

In the present CPI system, the shear or the separation of the two considered light rays d is assumed to be very 
small in comparison with the dimensions of the field of view. The right-hand side of equation (1) can be developed 
into a Taylor series around n(xy,z)- Ignoring the higher order terms and taking into account only the linear term 
(because d is small) yields 

Z2 a 

A l(x,y) = dj —n(x,y,z)dz (2) 

Zl 


The quantity A l(x,y) is measured in the recording plane in the form of bright interference fringes that appear where 


— Al(x, y) = 0, ±1, ±2 , ... 


(3) 
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with X being the light wavelength. The equation for the dark fringes is 




(4) 


Note that the measured signal, the left-hand side of equations (3) and (4), is two dimensional whereas the quantity to 
be determined n(x,y,z) is three dimensional in the general case. A direct measurement of n is possible for a plane 
refractive index field n(x,y). 

The dependence of the measured signal on the gradient of the refractive index, as in the techniques that rely on 


the measurement of the deflection angle € 


n r -2 1 dn 

~ J ~^ dz (described in ref. 6), can be seen in equation (2), where 


n a is the refractive index of air. The measurement of dn/dx is also possible by rotating the optical system 90° around 
the optical axis (--axis in this case). From equation (2) it follows that an interferogram with dn/dy will appear uni- 
formly illuminated; all light waves arrive in the recording plane with the same optical phase. This case is called the 
“infinite fringe mode’’ alignment of the interferometer. It can be aligned in a different way so that a set of equidistant 
parallel interference fringes appear for the object with dn/dy = a constant. 

For the study reported herein, the interferometer was set in the finite fringe mode, and the prism was rotated 
around the optical axis z to create a small angle between the fringes and the diffusing interface under investigation, 
making it easy to illustrate the concept that the gradient of the index of refraction is directly proportional to the con- 
centration gradient in the fluids (ref. 8). The final relations are given in the following equations (ref. 8). 

For n - n(y). 


And similarly for n = n(x). 


dn _ ^ dc 
dy dy 


dn 

dx 



(5) 


(6) 


where parameters K and Q are constant values that are proportional to the wavelength and inversely proportional to 
the object length L, the separation distance d , and the angle between the interface of miscible fluids under measure- 
ment and the fringe direction. Equations (5) and (6) do show that the concentration gradient is proportional to the 
refractive index gradient. By using this concept and assuming that the diffusion process is one dimensional and is 
characterized by a constant value of the diffusivity as described in the section Diffiisivity Calculation, the diffusivity 
for two miscible liquids in contact was obtained. 


Test Apparatus and Procedure 

Upon using the concept expressed in equations (5) and (6) from reference 8, the CPI measures the gradient of 
the index of refraction in a certain direction within a test cell. A Wollaston prism was used to image the sheared 
light beam by the concentration gradient near a miscible interface. A collimated, polarized beam of light from a low- 
power laser is passed through a transparent test cell. The light is then focused on the Wollaston prism, which splits 
the light into two beams that are slightly displaced from each other. When the beams are recombined, they produce 
interference fringes that indicate gradients of refraction in the test medium. The interferometer is arranged in its 
finite fringe mode where equidistant, parallel interference fringes appear in a uniform index-of-refraction field in the 
test cell. This is the case when only one of the fluids is present in the test cell. When a second liquid is introduced 
and diffusion occurs, the deviation or shift of a fringe from its undisturbed location is a measure of the index-of- 
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Figure 1. — Test apparatus, (a) Common-path interferometer (CPI), (b) Expanded shearing unit. 


refraction gradient within the test cell. It has been assumed that the index of refraction within the test cell is only a 
function of the vertical coordinate. The light beam is passed through the test cell horizontally, and it has been further 
assumed that each ray of light traverses a path of constant refractive index without refraction. Figure 1 shows a sche- 
matic of the CPI and the interferometer unit. The test cell used for the diffusivity measurement of a pair of miscible 
liquids is a quartz container 10 by 10 mm in cross section and 45 mm in height (fig. 2). A charge-coupled-device 
(CCD) camera captured the interferograms of the CPI and a time-lapse super video recorder (S-VHS) stored them 
on a tape for later analysis. Figure 3 shows a typical fringe image near the diffusing interface after the fluids came in 
contact. Figure 4 illustrates the trace of a typical fringe near the interface under investigation and also shows the data 
points that are easily identified and used for calculating the diffusivity of the liquids under measurement. These 
points are selected because they simplify the diffusivity calculations when an average of the diffusivity over some 
length of time is needed (explained in the next section). The fluids used in this measurement are silicone oils of vari- 
ous viscosities (table I). 


Diffusivity Calculation 

When two miscible liquids are put in contact, the two fluids are assumed to be pure and therefore only binary 
diffusion occurs. Let pjC*,/) and p 2 (x,t) represent the mass concentration of the two fluids at a vertical location * and 
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Figure 3. — Typical fringe near interface 30 min after 
liquids came in contact. (Note: The rest of fringes 
near the interface are not shown because their 
spacing is larger than the shift peak at this moment.) 


at a given time t. Far away from the interface, p loo and p 2oo are the densities of the two fluids. It is convenient to 
define the dimensionless molar concentration of liquid 1 as c — where M, is the molecular weight of 

MlPloo 

liquid 1. A similar dimensionless concentration for liquid 2 may be defined. For stationary liquids, the governing 
equation for diffusion is (ref. 9) 

^ = D^r (7) 

dt dx" 


Strictly speaking, this equation is valid only when the binary diffusion coefficient D and the local density in the liq- 
uid are constant. The coefficient of diffusion from our experiment will be extracted assuming that equation (7) is 
valid. The boundary conditions are 
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c — > 1 


as x — > oo 


and c — > 0 as x — > -©o 


(8) 


The solution for the concentration gradient is 


3c - 1 (9 

dx'ljKDt ( 

The dimensionless concentration gradient profile for liquid 2 is identical to that for liquid 1 . Assuming that the 
refractive index is linearly proportional to the concentration of the two fluids, the gradient of the refractive index has 
a similar profile. 


dc 

dx 



Figure 4. — Typical trace of fringe near interface after liquids came in contact. 
Data points for diffusivity D calculation are also shown on traced shifted 
fringe. (Note: All dimensions are normalized by the value of concentration 
gradient m.) 


TABLE I. — COMPARISON OF AVERAGED DIFFUSIVITIES 


Silicone oil liquid pairs, 
cSt 

Averaged diffusivities, cm 2 /s 

Technique 

T>l-(l/c) 

£>l/e 

•D Wiener 

1000/100 

(4.36±0.30)xl0 -8 
(5.1 l±0.73)xl0 -8 

(4.79±0.22)xl0 -8 

(5.41±0.75)xl0 -8 

(1.69±0.06)xl0 -8 

Wiener (ref. 10) 
Common-path interferometer 

1000/10 

(2.61±0.10)xl0- 7 

(2.83±0.72)xl(r 7 

(2.47±0.10)xl0- 7 
(2.34±0.68)x 1 0' 7 

(3.14±0.15)xl0- 7 

Wiener (ref. 10) 
Common-path interferometer 

1000/1 

(4.17±0.23)xl0“ 6 

(2.65±0.32)xl0~ 6 

(2.48±0.21)xl0 -6 

(1.93±0.29)xl0 -6 

(1.93±0.08)xl0 -6 

Wiener (ref. 10) 
Common-path interferometer 
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The maximum value of the concentration gradient is at x = 0 and is given by 


m = 



1 


( 10 ) 


For positive x, let x ; = — <5 ( {/ \ at the location where — = ( 1 - — V 7 . The solution given for the 
2 l ~\/e) dx V ey 

gradient can then be used to obtain the diffusivity: 


concentration 



(ID 


Similarly, let x = — by at the location where — = — m. Then D 1/e is obtained to be 

2 /e 0* 


T 

Dy = -A 
* I6t 


( 12 ) 


The values of 5 1 _ (1/e) and 8 1/e have been obtained from the measured refractive index gradient profiles. The time t 
after the fluids came in contact and each value of Sj_^ 1/e) or 8 1/e extracted from the recorded fringes (fig. 4 shows a 
typical fringe profile with the data points of interest on it) can be used to calculate the diffusivity from equations 
(11) and (12), respectively. For an average diffusivity, the values of 8 1-(1/e) and 5 1/e for various values of t were 
used and 8 2 was plotted versus t. The values for diffusivity, with its 95-percent confidence interval, were obtained 
by using linear regression and finding the best fitting straight line to such plots. 


RESULTS 

Table I shows the measured diffusivity (averaged over time) for several silicone oil liquid pairs. For compari- 
son, it also contains a set of data similar to that obtained with other techniques (ref. 10). The diffusivity was obtained 
with the assumption that the simple theory of binary diffusion holds. In all pairs of liquids tested, from the regres- 
sion analysis, a 95-percent confidence interval for the measured diffusivity is reported herein. The diffusivity meas- 
ured by the original Wiener’s technique reported in reference 10 appears to be the least accurate, perhaps because 
the refractive index must be known more precisely for the wavelength of the light used. In the case of diffusion 
across the interface of one silicone oil pair (1 and 1000 centistokes), the difference in diffusivity obtained by the 
various techniques is large (36 percent because of the nongaussian distribution in the concentration gradient across 
the interface). This pair has the largest density difference (Ap/p = 0.19) among the various techniques compared 
herein. One possible explanation for the gap between data sets is that the simple model of diffusion that was 
assumed may not hold for this pair. 

The reliability of the CPI has made it attractive for application in diffusion measurements near moving inter- 
faces of miscible fluids in selected space experiments (ref. 11 ). 
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CONCLUSIONS 


A common-path interferometer (CPI) was developed and used in its finite fringe mode. A novel approach was 
taken to provide a fringe profile that resembled the solution to a simple model of a diffusion equation. The approach 
involved setting a small angle between the fringe orientation and the interface of the two miscible liquids. Then a 
mathematical model of the optics for this situation was given. From the analysis of the fringe profile obtained in this 
manner, the diffusion coefficients (diffusivity) of pairs of miscible fluids were calculated. 

The technique used herein was compared with Wiener’s original technique, which requires the values for the 
refractive indices of the two miscible fluids. His appears to be the least accurate because the refractive indices must 
be known more precisely for the wavelength of light used. This information is not required for the CPI technique. 
The diffusion coefficients, measured by interferometry and a modified Wiener’s technique that analyzes the 
measured concentration gradient profiles, are in good agreement. The exception is the case of diffusion across the 
interface of a silicone oil pair of 1 and 1000 centistokes, for which the difference is as much as 36 percent. This pair 
has the largest density difference. An explanation might be that the simple model of diffusion used herein to extract 
the diffusion coefficient did not strictly hold for this pair. The measured concentration gradient profile is not a 
gaussian symmetric profile. Also the diffusivity may be dependent on the concentration itself, which requires further 
analysis of the data to extract such information from similar interferometric signals. 

The approach reported herein and the application of the technique to double-diffusive problems will be the sub- 
ject of future studies. The present CPI technique has been adopted for diffusion measurement near moving interfaces 
of miscible fluids in a selected space experiment in the near future. 
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